BAR4R-CONF-04/04 
SLAC-PUB-10593 
August 2004 



Measurements of yl+ Branching Fractions of Cabibbo-Suppressed 

Decay Modes. 

The BABAR Collaboration 
February 7, 2008 



Abstract 

We have measured the branching fractions of the Cabibbo-suppressed decays — > A^K^ and 
A'^ Z!^K~^ relative to the Cabibbo-allowed decay modes A'^ yl'^vr"^ and Af U^7t~^ to be 
0.044 ± 0.004 (stat.) ± 0.002 (syst.) and 0.040 ± 0.005 (stat.) ± 0.004 (syst.), respectively We 
also present the first observation of A'^ A^K^tt^tt^ and have measured the branching fraction 
relative to A+ A^tt+ to be 0.266 ± 0.027 (stat.) ± 0.032 (syst.). The upper limit of the 
branching fraction into the decay A'^ — > E^K^tt^tt^ relative to A^ — > U^tt~^ has been measured 
to be < 3.9 X 10~^ at the 90% confidence level. This analysis was performed using a data sample 
of 125 fb~^ (integrated luminosity) collected by the BABAR detector at the PEP-II asymmetric- 
energy B Factory at the Stanford Linear Accelerator Center. All results presented in this conference 
contribution are preliminary. 
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1 INTRODUCTION 



Beginning with the first observation of the charmed baryon A'^ in 1979 by MARK-II and BNL 
^ our knowledge of the physics of charmed baryons developed less rapidly than that of the 
charmed mesons. This is due to the smaher baryon production cross section, shorter life time, 
and, in e^e~ storage rings, the absence of a cleanly observable Ac Ac resonance. During the last 
few years there has been significant progress in the experimental study of the hadronic decays 
of charmed baryons. Recent results on masses, widths, lifetimes, production rates and the decay 
asymmetry parameters have been published by different experiments; among them the discoveries 
of Cabibbo-suppressed decays A+ pcj) by CLEO 3 , and A+ A'^K^ , A+ S^K+ by Belle 

in- 

The precision in the measurements of branching fractions is only about 40% for many Cabibbo- 
favored modes while for Cabibbo-suppressed decays the precision is even worse. As a conse- 
quence, we are not yet able to distinguish between the decay rate predictions made by different 
models, e.g., the quark model approach to non-leptonic charm decays and the Heavy Quark Effec- 
tive Theory(HQET) [SICZIIHI- 

In this paper we present a study of A'^ baryons produced in the e^e~ — > qq continuum at 
BABAR. We present improved measurements of the Cabibbo-suppressed decays A'^ A^K^ and 
S^K^, report the first observation of A'^ A^K'^ir^n'^ , and set an upper limit on A^ 
vr+vr . Here and throughout this paper, inclusion of charge-conjugate states is implied. 

2 THE BABAR DETECTOR AND DATASET 

The data sample used in this analysis consists of 125 fb~^ integrated luminosity recorded between 
October 1999 and June 2003 with the BABAR detector at the SLAC PEP-II storage ring. The PEP- 
II facility operates nominally at the T(4S') resonance, providing collisions of 9.0 GeV electrons on 
3.1 GeV positrons. The data set includes 112 fb~^ collected at the T{AS) ( "on-resonance" ) and 13 
fb"^ collected below the BB threshold ("off- resonance i . 

A detailed description of the BABAR detector can be found elsewhere [0] ; only detector compo- 
nents most relevant to this analysis are mentioned here. Charged-particle trajectories are measured 
by a five-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH), oper- 
ating in the field of a 1.5-T solenoid. Charged particles are identified by combining measurements 
of ionization energy loss {dE/dx) in the DCH and SVT with angular information from a detector 
of internally refiected Cherenkov light (DIRC). Photons are identified as isolated electromagnetic 
showers in a CsI(Tl) electromagnetic calorimeter. 

We have used Monte Carlo simulations (MC) of the BABM detector based on GEANT4 (THl 
to optimize our selection criteria and to determine signal efficiencies. These simulations take into 
account the varying detector conditions and beam backgrounds during the data-taking period. 

3 ANALYSIS METHOD 

For this analysis, the particle identification is important. A track is identified as a kaon, pion, or 
proton if it is projected to pass through the fiducial volume of the DIRC and the reconstructed 
cone of Cerenkov light is consistent in time and angle with the measured track momentum. This 
information is augmented with dE/dx measured with the SVT and DCH. Photons are detected in 
the Csl calorimeter. 
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Candidates for are reconstructed in the decay mode — > piT~ . We fitted the p and tt~ 
tracks to a common vertex and required the probabihty of of vertex fit to be greater than 
0.1 %. We also required the (three-dimensional) flight distance of each A^ candidate between its 
decay vertex and the primary vertex to be greater than 0.2 cm. We fitted the invariant mass 
of A^ candidates, with a double-Gaussian function of common mean to represent the signal and 
a second-order polynomial to represent the background. This fit is shown in Fig. ^ The fitted 
resolution is a = cjrms = 1-5 MeV/c^, where cjrms is defined by 

2 — r 2 : r 2 
C^RMS = /I'^l + /20-2> 

where /i and /2 are fractions of signal yield corresponding to Gaussian functions one and two, re- 
spectively, and (Ti and (T2 are the two corresponding widths. We required the mass of A^ candidates 
to be in the range 1113 MeV/c^ < M^o < 1119 MeV/c^. The candidates were reconstructed 
in the decay mode — > A^'y using the A^ sample and photons with a calorimeter cluster energy 
greater than 0.1 GeV. The mass difference (Myjo^ — Mj^o) is shown in Fig. [21 Fitting with two 
Gaussian functions of common mean for the signal contribution, and with a third-order polyno- 
mial for background, we obtained a resolution a = ctrms = 4.0 MeV/c^ and a mean ( difference 
between the invariant masses) of 767 it 1 MeV/c^. We accept candidates with My^o^ — Mj\o within 
± 10.0 MeV/c^. (2.5 a) of the measured mean value. 

To suppress combinatorial and BB backgrounds, we required A'^ candidates to have scaled 
momentum Xp = p*/Pmax > O-^i here p* is the reconstructed momentum of the A'^ candidate in 
the e"'"e~ center of mass, p^ax = a/s/4 — M^, y/s is the total center of mass energy and M is the 
reconstructed mass of the A'^ candidate. The signal detection efficiency is obtained from MC with 
particle identification efficiency corrections based on data. 




'«(p.)(GeV/c') 

Figure 1: The invariant mass oi p-K combinations (GeV/c^). 
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4 PHYSICS RESULTS 

4.1 Measurement of the decays yl^ 



A^K+ and At E^K^ 



The Cabibbo-suppressed decay A^ — > A^K'^ was first measured by the Belle Collaboration j^. 

For our analysis, the A^ and were combined to form as described in Sec. 01 Invariant 
mass distribution of AS'K^ combinations is shown in Fig. |3J The mass distribution was fitted with 
a Gaussian function for the signal, and a second-order polynomial for combinatorial background. 
We obtain a raw yield of 1164 it 107 (stat.) events, with a 10.9 standard deviations significance 
and with fitted width a = 5.5 it 0.7 MeV/c^, consistent with the MC prediction of 6.0 MeV/c^. 

For normalization, we used the decay A^ — > /I^tt"*". The invariant mass distribution of APtt'^ 
combinations is shown in Fig. At mass values below the peak centered at the A^ mass a broad 
distribution at 2.2 GeV/c^ is visible which is a refiection due to A^ — > S^'k'^ with a missing 
7. Additionally, at 2.3 GeV/c^ we see a shoulder, the upper edge of a S'c reflection whose full 
shape extends through the entire A^ signal region. The distribution was fitted using two Gaussian 
functions with same mean for signal, a square wave function for each reflection, and a 7th-order 
polynomial for combinatorial background. We obtained the width a = 8.2 MeV/c^, which is 
consistent with the experimental resolution of about 8.0 MeV/c^, and a raw yield of 33594 it 367 
(stat.) events. We calculate the ratio B{A't^ K^) /B{At ^ yl'^vr^) in on-resonance and off- 
resonance data. The results are comparable within uncertainties, so off-resonance data is combined 
with on-resonance data. Using signal MC, the relative signal reconstruction efficiency is found to 
be e(yl+ ^ A^K+) / e(yl+ ^ A^tt+) = 0.781 ± 0.004 (stat.). With this value we calculate: 



B{A+ A^K- 
B{A+ ^ #7r+ 



0.044 ± 0.004 (stat.) ± 0.002 (syst.) 



We provide a detailed description of the sources of systematic uncertainty for this and other mea- 
sured decay modes in Sec. EJ 

The Cabibbo-suppressed decay A^ —>■ S^K~^ was first measured by the Belle collaboration [1], 
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Figure 3: The invariant mass of APK^ combinations (GeV/c^) 
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Figure 4: The invariant mass of APit~^ combinations {GeV /(?] 



under the restriction of the scaled momentum to the region Xp > 0.6. For our analysis, we com- 
bined a and a candidate to form , requiring Xp > 0.5 as before. The invariant mass of 
E^K^ combinations is shown in Fig. El We improved the invariant mass resolution by plotting the 
corrected mass e.g., M^ok+ — Mjjo + M^P^ , instead of M^oj^+ where M^o is the reconstructed 
mass of and M^P^ is the mass of from PDG We fit the distribution using a Gaussian 
function with a fixed width a = 6.0 MeV/c^ (as determined from MC simulations) for signal and a 
third-order polynomial for combinatorial backgrounds. The fit yields 387 it 48 (stat.) events with 
a 8.1(7 statistical significance for A'^ baryons having decayed to U^K^. For normalization, we use 
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the Cabibbo-allowed decay mode Af — > E^n'^ . The invariant mass of U^ir^ combinations is shown 
in Fig. ini The fit uses a Gaussian function for the signal and a third-order polynomial for back- 
ground. The measured width a = 6.7 it 0.1 MeV/c^, which is consistent with the MC prediction 
of cr = 7.0 MeV/c^. The fitted yield is 12450 it 170 (stat.) events. The relative reconstruction 
efficiency is measured to be e{A+ U'^K+) / e{A+ U°tt+) = 0.780 ± 0.001 (stat.) using 
signal MC. The resulting relative branching ratio is 



mi 



0.040 ± 0.005 (stat.) ± 0.004 (syst.) 
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Figure 5: The invariant mass of S^K~^ combinations (GeV/c^). 



4.2 Measurement of the decay /l^ A'^K'^n^n 

To measure the Cabibbo-suppressed decay A'^ — > A^ K^tt^tt^ we used the selection criteria 
described in Sec. 01 but with scaled momentum restricted to Xp > 0.6. The invariant mass distri- 
bution of A^K~^TT^TT~ combinations is shown in Fig. [Tj The distribution is fitted with a Gaussian 
function for the signal shape and a third-order polynomial for background. We obtained the width 
0- = 9.9 ± 1.0 MeV/c^ and signal yield of 2591 ± 258 (stat.) events for the A+ A'^K+tt+tt- 
decay with lO.lcr statistical significance. A small fluctuation is seen at the mass 2.26 GeV/c^. 
The effect of this fluctuation on the signal yield has been included in the fitting systematic. For 
normalization mode we use A'^ — > A^7r~^ with scaled momentum at Xp > 0.6, for which we obtained 
a raw yield of 22173 ± 287 (stat.) events. The relative signal reconstruction efficiency is measured 
to be e(yl+ A'^K+tt+tt-) / e(yl+ yl°7r+) = 0.442 ± 0.004 (stat.). The resulting branching 
ratio is 

^ , , ;n— n = 0.266 ± 0.027 (stat.) ± 0.032 (syst.). 

This is the first measurement of the Cabibbo-suppressed decay yl+ — > A^K'^tt'^tt" . 
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Figure 7: The invariant mass of A^K'^tt'^t: combinations (GeV/c^) 



4.3 Search for the decay of A+ S^K+n+n- 

We searched for the decay A'^ E^K^tt^-k^ using the selection described in Sec. El and 
restricting the scaled momentum to Xp > 0.6. The fit of the invariant mass distribution (Fig. [H|) 
yields 41 ± 51 (stat.) events for the A+ E^K+ TT'^TT decay with 0.8a statistical significance. The 
width a = 10.0 MeV/c^ and mean = 2285.0 MeV/c were fixed to values obtained from MC. Using 
the decay mode A'^ E^7r~^ for normalization, we find a raw yield of 8785 it 131 (stat.) events 
for the decay A^ S^Tr~^ at Xp > 0.6. The relative reconstruction efficiency is measured to be 
e{A+ E^K+TT+TT-) I e(yl+ i:°7r+) = 0.390 ± 0.002 (stat.). We do not observe any significant 
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signal for U^K^tt'^tt . Therefore, we calculate the upper limit using the Feldman and 

Cousins method including systematic uncertainties. We find: 

B{A+ U'^K+TT+TT-' 



< 3.9 X 10" 



vr ' 



CL. 



5 SYSTEMATIC STUDIES 

We have considered several possible sources of systematic uncertainties in our measurements. The 
systematic uncertainty due to each cut used in candidate selection is ~ 1 %. The photon spectrum is 
different in measured and reference decay modes, so we studied the effect by changing the photon 
energy and assigned a systematic uncertainty of 2.8 % to the branching ratio {A'^ — > X!^K~^ / 
At - U'7r+). 

We have studied possible biases due to our fitting procedure. The shape of background has 
been varied by changing the degree of the polynomial function as well as the shapes of U and 
Sc reflections (in case of A^ — > A^tt~^ decay). Any change in signal yields is being taken as a 
systematic uncertainty. We also varied the signal width (cj) by one standard deviation, with an 
change in yields interpreted as a systematic uncertainty. We assign a systematic uncertainty due to 
fit bias of 5.1 % for A+ A^K+ / A+ ^ A^tt+ , 8.0 % for A+ E^K+ / A+ ^ S^tt+ , and 10.1 % 
for At — > A^K^'k^t:~ / A^ AS'tt'^ . The systematic uncertainty associated with the fitting is the 
dominant one. 



6 SUMMARY 

We report on a measurement of the branching ratio of the Cabibbo-suppressed decays A^ A^K~^ 
and At U^K^ with improved accuracy. We also report the first observation of the Cabibbo- 
suppressed decay At — > A^K'^ir^ir' , and we set an upper limit on the At — > U^K^tt^tt~ decay. 
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The results for these decay modes are summarized in Tabled All results reported in this paper are 
preliminary. The expectations from the quark model ^6 are B{A'^ — > A^K^)/B{A'^ — >■ yl^vr"'") = 
[0.039-0.056] and B{A+ U'^K+)/B{A+ i;°7r+) = [0.033-0.036]. The results are in agreement 
with the predictions. 
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Signal Mode 


Signal 
Yield 


Reference 
Mode 


Reference 
Yield 


Relative 
Efficiency 


'-^signal 

reference 


Other Measurements 




1164 ± 107 


AOtt+ 


33594 ± 367 


0.781 ± 0.004 


0.044 ± 0.004 ± 0.002 
{xp > 0.5) 


0.074 ±0.010 ± 0.012 4, 
(xp > 0.5) 




2591 ± 258 


A^n+ 


22173 ± 287 


0.442 ± 0.004 


0.266 ± 0.027 ± 0.032 
(xp > 0.6) 






387 ± 48 




12450 ± 170 


0.780 ± 0.001 


0.040 ± 0.005 ± 0.004 
{xp > 0.5) 


0.056 ± 0.014 ± 0.008 4 
{xp > 0.6) 




41 ±51 


U^7r+ 


8785 ± 131 


0.390 ± 0.002 


< 3.9 X 10-2 @ 90% CL 
{xp > 0.6) 





Table 1: Summary of results obtained in this analysis. The last column shows the current other measurements of each decay mode, 
where available. 



